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ABSTRACT 

The bolometric luminosities and black hole (BH) masses are estimated by various methods for a sample of 
young radio galaxies with known ages. We find that the ages are positively correlated with the bolometric lumi- 
nosities in these young radio galaxies, which is consistent with theoretical prediction based on radiation pressure 
instability of accretion disk in Czerny et al. The ages of young radio galaxies are also found to be consistent with 
the theoretical durations of outbursts in BH mass and accretion rate (in Eddington unit) plane, where the outbursts 
are assumed to be triggered by the radiation pressure instabilities. Our results provide the observational evidence 
for the radiation pressure instability, which causes limit-cycle behavior, as a physical mechanism that may be 
responsible for these short-lived young radio galaxies. 

Subject headings: galaxies: active, quasars — accretion, accretion disks — black hole physics, instabilities 
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1. INTRODUCTION 

Gigahertz-peaked spectrum (GPS: projected linear size D < 
1 kpc) and compact steep-spectrum (CSS: projected linear size 
D < 20 kpc) radio sources constitute a large fraction (~ 40%) 
of the powerful radio source population. Their radio spectra are 
simple and convex with peaks close to 1 GHz and 100 MHz 
for GPS and CSS sources respectively (see O'Dea 1998 for a 
review). Two basic models have been proposed to explain the 
CSS/GPS phenomena: (1) youth scenario: the radio sources 
are compact because they are young, and the jets will expand 
into large-scale 'classical' lobes during the course of their life 
(e.g., Fanti et al. 1995); (2) frustration scenario: the radio 
emitting plasma is permanently confined to the region within 
the host galaxy by a cocoon of dense gas and dust for their 
entire life (e.g., van Breugel 1984). Both kinematic ages esti- 
mated from proper motions of hot spots (e.g., Giroletti & Po- 
latidis 2009 and references therein) and radiative ages deduced 
by synchrotron theory based on the radio spectra of lobes (e.g., 
Murgia 1999 and reference therein) provide values from ~ 10 2 
to ~ 10 5 years, which strongly support the youth scenario. 

Duty cycle is one of key problems in understanding the evo- 
lution of active galactic nuclei (AGN, e.g., Wang et al. 2009, 
and references therein). The extended radio emission in galax- 
ies and quasars provides us with an opportunity to probe their 
history through the ages derived from the structural and spec- 
tral information of their radio lobes. Both the young compact 
symmetric object (CSO) with a much older extended relic emis- 
sion (e.g., Owsianik et al. 1998) and the 'double-double' radio 
source which consists of a new lobe-pair near the core and a 
more distant faded lobe-pair (e.g., Lara et al. 1999; Schoen- 
makers et al. 2000) can be explained in terms of recurrent activ- 
ity in their nuclei. These possibly re-born radio sources may be 
intermittent on timescales of ~ 10 4 - 10 5 years (e.g., Reynolds 
& Begelman 1997). 

A standard, geometrically thin accretion disk (SSD, Shakura 
& Sunyaev 1973) is known to be unstable in the innermost, ra- 
diation pressure dominated regions when accretion rate is high 
(e.g., Lightman & Eardley 1974). The thermally unstable in- 
ner region will possibly lead to the so-called limit-cycle be- 
havior, which has been confirmed by numerical simulations 
of the time-dependent disk (e.g., Honma et al. 1991; Li et 



al. 2007). On the observational side, a perfect candidate of 
black hole (BH) X-ray binary GRS 1915+105 has been known 
to show the theoretically predicted limit-cyclic luminosity vari- 
ations (e.g., Belloni et al. 1997; Nayakshin et al. 2000; Janiuk 
et al. 2002; Merloni & Nayakshin 2006). Czerny et al. (2009) 
found that the outburst durations are generally from ~ 10 2 to 
~ 10 s years for BH mass M B h = 10 7 - 10 10 M Q and viscosity 
parameter a = 0.02-0.2, which are roughly consistent with the 
observed ages of young radio galaxies, where the viscosity law 
with a geometrical mean between the gas and the total pressure 
for scaling of the stress, i.e., T r( f, oc 0P ga sAot) 1//2 , is adopted when 
modeling the time evolution of SSD under the radiation pres- 
sure instability. An important empirical correlation between the 
outburst duration, 7b urst , bolometric luminosity, Lboi, and vis- 
cosity parameter, a, is given by Czerny et al. (2009), 
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Wu (2009) found that most of young radio galaxies have 
relatively high accretion rates with average Eddington ratio 
< logLboi/^Edd >= -0.56 for a sample of 51 sources (32 CSS 
+ 19 GPS), where L Edd = 1.38 x 1O 38 M BH /M is Eddington 
luminosity. The accretion disk in these high Eddington ratio 
CSS/GPS sources may have suffered the radiation pressure in- 
stability, which provide a possibility to explain their short-lived 
activities. In this Letter, we estimate the bolometric luminosi- 
ties and BH masses for a sample of young radio galaxies with 
known ages. We then test whether the ages of young radio 
galaxies are consistent with the durations of outbursts or not, 
where the outbursts are assumed to be triggered by the radia- 
tion pressure instabilities of accretion disk (e.g., Czerny et al. 
2009). Throughout this paper, we assume the following cos- 
mology: H = 70 km s _1 Mpc _1 , fi = 0.3 and A = 0.7. 

2. THE SAMPLE 

We compile a sample of 37 young radio sources from the data 
available in the literature. The kinematic ages are mainly se- 
lected from Giroletti & Polatidis (2009, and references therein), 
while the ages estimated from synchrotron cooling time are 
mainly taken from Murgia et al. (1999). We give the lower 
limits of kinematic age for four sources (0134+329, 0538+498, 
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1458+718, 2249+185) assuming the upper limit for the sepa- 
ration of hot spots to be the light speed c in this work, since 
their radiative ages may be much shorter than the actual source 
ages in these jets or hot spots dominated sources (see Murgia 
et al. 1999 for more details). The kinematic age of CSS source 
1328+307 (3C 286) is estimated in the present work from its 
projected liner size and jet propagation speed ~ 0.1c (An et al. 
2004). Several other individual sources with known ages are 
also selected from the literature. In total, our sample consists of 
14 GPS and 23 CSS sources (see Table 1 for more details). 

3. BOLOMETRIC LUMINOSITIES AND BH MASSES 

The bolometric luminosity of AGN is sometimes approx- 
imated from the optical luminosity, since the integration of 
spectral energy distribution (SED) is usually hampered by the 
lack of multi-wavelength observational data that spans many 
decades in wavelength. It is found that the bolometric lumi- 
nosity is roughly 9 times the optical luminosity, Lboi = 9Lsioo, 
where L5100 = AL A (5 100) (Kaspi et al. 2000). It should be noted 
that the monochromatic optical luminosity cannot be used to 
derive the bolometric luminosity if the optical luminosity is 
contaminated by the jet emission in radio loud AGN (e.g., BL 
Lacs) or obscured by the putative torus. The broad line emis- 
sion is believed to be illuminated by the ionizing luminosity 
from the accretion disk, and, therefore, it can be used as a good 
indicator of thermal optical emission. We estimate the bolomet- 
ric luminosities for 9 sources from their broad line luminosi- 
ties (e.g., H/3, C IV) using the empirical correlations between 
the luminosities of different emission lines and bolometric lu- 
minosities derived from a sample of radio-quiet AGN in Wu 
(2009). The bolometric luminosities of 2 sources fitted from 
their SEDs are selected directly (2342+821, 1458+718, Woo & 
Urry 2002). The bolometric luminosity of 1 source have been 
corrected for beaming effect is also included (1413+135, Xie et 
al. 2006) . The bolometric luminosities of 4 optically Type 1 
CSS sources are estimated from their optical luminosities with 
Lboi = 9L5100, since that, normally, the steep-spectrum sources 
do not suffer strong contamination from the Doppler-boosted 
emission of the jet due to the possible large inclination angles. 

In the standard model, narrow emission lines (e.g., [O III]) 
are also believed to be excited through photoionization by the 
UV continuum produced by the accretion disk (e.g., in radio- 
quiet AGN). However, it is more complicated in young radio 
galaxies, since shocks caused by the interaction between the 
jet and interstellar medium (ISM) may also be important for 
the morphology, kinematics, and the excitation of narrow emis- 
sion lines (e.g., Labiano et al. 2005; Holt et al. 2006). Moy 
& Rocca-Volmerange (2002) proposed that the shock ioniza- 
tion may be only important for sources with intermediate ra- 
dio sizes (2 kpc<Z)<150 kpc, e.g., CSS sources), while AGN 
photoionization is dominant in the most compact sources (D<2 
kpc, e.g., GPS sources) and the extended sources with D>150 
kpc. Moreover, the jet sizes in all GPS sources of our sample 
are less than 0.2 kpc (see Table 1), which are much less than 
the kpc-scale of narrow line regions. This also support that the 
shock contribution to narrow lines from the jet-cloud interac- 
tion may be not important and, therefore, narrow line emission 
may still be a good indicator of central thermal emission in 
these GPS sources. To derive the bolometric luminosities for 
more sources, we investigate the relation between [O III] lumi- 
nosities and bolometric luminosities for a sample of radio-quiet 
Type 1 AGN, which include 20 Seyferts and 23 low red-shift 



QSOs. The Seyferts are selected from Kaspi et al. (2000), 
where the [O III] luminosities are chosen from Kraemer et al. 
(2004), and the QSOs are selected from Alonso-Herrero et al. 
(1997). To be consistent with other methods, we calculate the 
bolometric luminosities from Lboi = 9L5100 for these radio-quiet 
AGN. The relation between the bolometric luminosities and [O 
III] luminosities is shown in Figure 1, and its linear regression 
is, 

logL bol = 0.95 ± 0.081ogL [O mj + 5.39 ±3.17. (2) 

We find that bolometric correction factors, Lboi/L[0 in], are 
about 2000 and 3400 for 20 Seyferts and 23 QSOs respectively. 
The bolometric luminosities for 12 GPS and 9 CSS sources 
were calculated with Eq. (2) from their [O III] luminosities. 
One should keep in mind that the bolometric luminosities of 9 
CSS sources may be over-estimated due to the possible shock 
contribution to [O III] luminosities (see discussions in Sect. 5). 
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FIG. 1. — The correlation between bolometric luminosities and [O III] lumi- 
nosities for a sample of radio-quiet Type 1 AGN. The circles and squares are 
correspond to 20 Seyferts and 23 QSOs respectively. The solid line is the best 
fitting. 

Wu (2009) have estimated BH masses for a sample of 65 
young radio galaxies from different methods. We estimate BH 
masses of other 14 young radio galaxies, which are not included 
in Wu (2009), from the relation between host galaxy absolute 
magnitude M R in /?-band and BH mass proposed by McLure et 
al. (2004), 

log 10 (M B H/M Q ) = -0.5M s -2.74. (3) 

4. RESULTS 

We estimate the bolometric luminosities for 37 young radio 
sources by various methods. The relation between ages and 
bolometric luminosities for young radio galaxies is presented 
in Figure 2, where the bolometric luminosities of CSS sources 
calculated from [O III] luminosities are marked with extra open 
squares. We find that the ages are positively correlated with the 
bolometric luminosities in these young radio galaxies, which is 
consistent with the prediction of radiation pressure instability 
in accretion disk (Czerny et al. 2009). We plot the dashed lines 
corresponding to Eq. (1) with viscosity parameters a =10" 4 and 
1 in Figure 2. We find that most sources are locate to the left 
of the line with a = 1 (Lboi-Tbm-st), where the a = 1 line gives 
an upper limit of the theoretical outburst duration. The relation 
between bolometric luminosity and half of the outburst dura- 
tion with a = 0.1 is also presented in Figure 2 (Lboi-Tburst/2, 
solid line), which is converted from the Eq. (1). We find that 
the real ages of sources in our sample are roughly consistent 
with the theoretical prediction of 7b urst /2 with a ~ 0.1, since 
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that sources should have the real ages between and 7b ulst (e.g., 
0<age<7b urst ), and, on the average, the real ages should corre- 
spond to rb ulst /2 with some dispersion around this value. 
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FIG. 2. — The correlation between ages and bolometric luminosities for 
young radio galaxies. The dashed lines are theoretical predictions from the 
disk instability model for Lbol - Thurst with a =10~ 4 (upper) and 1 (lower) re- 
spectively, while the solid line is for — T^ urst /2 with a =0 . 1 . The bolometric 
luminosities of CSS sources calculated from [O III] luminosities are marked 
with extra open squares. 
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FIG. 3. — Comparison of the ages for young radio galaxies with the theo- 
retical durations of outbursts triggered by disk instability in Mbh _ "i plane. 
The triangles, squares, pentagons, and circles correspond to the age in bins of 
[0-5 x 10 2 ], [5 x l() 2 -5 x 10 3 ], [5 x 10 3 -5 x 10 4 ], and [5 x 10 4 -5 x 10 5 ]yrs 
respectively. The solid lines are the half of theoretical outburst durations pre- 
dicted by the disk instability model with a = 0. 1 , which are converted from the 
case of a = 0.2 in Czerny et al. (2009, Fig. 5) using the scaling relation of Eq. 
(1). 

The BH masses are calculated from the luminosities of host 
galaxy, or collected from the literature. The Eddington ratios, 
m =Lboi/^Edd, were computed from the estimated bolometric 
luminosities and BH masses. We then compare the ages of 
young radio sources with the half of outburst durations of in BH 



mass-accretion rate (in Eddington unit, Mbh _ tn) plane, where 
the outbursts are assumed to be triggered by the radiation pres- 
sure instabilities of accretion disk. To do so, we convert the 
contour map of a = 0.2 in Czerny et al. (2009, Fig. 5) to the 
case of 7b urst /2 with a = 0.1 using the scaling relation of Eq. 
(1), since we find that most of the young radio sources can, on 
the average, be described by a ~ 0.1 in the Lbd-age relation. 
The ages of young radio sources are divided into four bins as 
[0-5 x 10 2 ], [5 x 10 2 -5 x 10 3 ], [5 x 10 3 -5 x 10 4 ], and [5 x 10 4 - 

5 x 10 5 ]yrs according to the theoretical lines of 10 2 , 10 3 , 10 4 
and 10 5 yrs of Czerny et al. (2009) respectively. The results are 
shown in Figure 3. We find that, on the average, the ages of 
young radio sources are also roughly consistent with the half 
of outburst durations with a ~ 0. 1 in Mbh _ in plane, even with 
substantial scatter. In particular, there are 12 sources in the bin 
of [5 x 10 3 -5 x 10 4 ]yrs (red-pentagons), which coincide with 
the theoretical prediction of 10 4 yrs (red-solid line) very well. 

5. SUMMARY AND DISCUSSION 

We estimate bolometric luminosities and BH masses for a 
sample of young radio galaxies with known ages. We find a 
positive correlation between the ages and bolometric luminosi- 
ties in these young radio sources (see Fig. 2). This positive 
correlation is roughly consistent with the theoretical prediction 
based on the radiation pressure instability of accretion disk in 
Czerny et al. (2009). We find that the real ages of most sources 
are shorter than the upper limit of outburst duration with a = 1 , 
and, on the average, correspond to the half of outburst durations 
with a ~ 0.1 (see Fig. 2). We further compare the ages with 
the half of outburst durations in M^-in plane, and find that the 
observational data is also roughly consistent with the theoreti- 
cal prediction with a ~ 0.1 (see Fig. 3). Therefore, our results 
provide observational evidence for the radiation pressure insta- 
bility as a physical mechanism that may also be responsible for 
these short-lived young radio galaxies, beside the perfect appli- 
cation to the brightest Galactic X-ray binary GRS 1915+105. 

The bolometric luminosities for young radio galaxies are es- 
timated from their SED fittings, optical luminosities, broad line 
luminosities, or narrow line luminosities. In these methods, the 
most controversial one is that derived from the [O III] lumi- 
nosity, due to the possible shock contribution caused by the 
jet-ISM interaction. However, the shock contribution may be 
not important in the GPS sources as suggested by the ioniza- 
tion diagnostics of different narrow emission lines (e.g., Moy 

6 Rocca-Volmerange 2002). Wu (2009) also found that ac- 
cretion activities may still play an important role in shaping 
the kinematics of [O III] narrow line in young radio galaxies. 
The bolometric luminosity calculated from [O III] luminosity is 
logL bol ([0 III]) = 45.12 for the GPS source 1404+286, which 
is consistent with that calculated from the broad line luminos- 
ity log Lbol (H/3) = 45.19 very well, where the [O III] luminosity 
is selected from Marziani et al. (2003). This also support that 
the shock contribution to [O III] luminosity may not be impor- 
tant, and the [O III] luminosity can still be regarded as a good 
indicator of central thermal emission in these small radio-sized 
GPS sources. Labiano et al. (2005) found that the contribu- 
tion to narrow line luminosity from the shocked gas is between 
30 and 70% for CSS source 3C 303.1, and this ratio is even 
lower in other two CSS sources (3C 67 and 3C 277.1), based 
on the ionization diagnostics of different narrow emission lines. 
Therefore, the bolometric luminosities of 9 CSS sources calcu- 
lated from the [O III] luminosities may not be over-estimated 
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too much, which will not change our main conclusion. 

The limit-cycle behavior triggered by the radiation pressure 
instability of accretion disk provide a good explanation to the 
luminosity variations of X-ray binary GRS 1915+105 (e.g., 
Nayakshin et al. 2000; Janiuk et al. 2002; Merloni & Nayak- 
shin 2006). Czerny et al. (2009) proposed that the young ra- 
dio galaxies may have the similar origin, and the unstable disk 
caused by radiation pressure may be the physical reason for 
their short-lived radio activities. We note that the unique limit- 
cycle variability only appears in the brightest X-ray binary GRS 
1915+105 when it radiates at high luminosity (e.g., Belloni et 
al. 1997). Wu (2009) found that most of the young radio galax- 
ies have relatively high Eddington ratios, and the accretion disk 
may also have suffered the radiation pressure instabilities in 
these young radio sources. We collect a sample of 37 young 
radio galaxies, and explore whether these young radio sources 
are consistent with the predictions of disk instability model in 
Czerny et al. (2009) or not. We test this in both the Lbd-age 
relation and the ages in Mbh - m plane. We find that the real 
ages of most sources are, indeed, less than the upper limit of 
the theoretical outburst durations with a = 1 (see Fig. 2). The 
viscosity parameter, a, should be larger than 10~ 4 , since that the 
real ages of most sources are longer than its predicted outburst 
durations (see Fig. 2). We find that the real ages of these young 
radio sources are roughly consistent with the model predictions 
with a ~ 0. 1, if we assume the real ages, on the average, should 
correspond to the half of outburst durations (see Fig. 2). The 
real ages of the young radio sources are also found to be consis- 
tent with the half of outburst durations with a ~ 0.1 mM^-m 



plane (see Fig. 3). It should be noted that there are still uncer- 
tainties in estimating BH masses, bolometric luminosities and 
ages in these young radio galaxies. Moreover, the simplified 
disk model was used to derive the Eq. (1), where corona and jet 
are neglected (see Czerny et al. 2009 for more details). With all 
these uncertainties in mind, it seems that the observational data 
of young radio galaxies are in good agreement with the theoret- 
ical predictions of disk instability model. We note that the disk 
instability model also predicts a complex behavior of the jet ac- 
tivities depending on the accretion rate and BH mass: (1) con- 
tinuous jet at low accretion rate where the accretion disk is sta- 
ble (e.g., m < 10~ 2 , low luminosity AGN); (2) intermittent com- 
pact jet with short outburst duration (< 10 4 yr) at high accretion 
rate, where the disk is unstable and the short-lived outburst of 
central activity cannot escape from host galaxy (e.g., m > 10~ 2 , 
GPS and CSO sources); (3) intermittent jet with long outburst 
duration (> 10 4 yr) at high accretion rate (e.g., m > 10~ 2 , CSS 
sources), where the disk is unstable and the jet is able to es- 
cape from the host galaxy and possibly grow to large-scale ra- 
dio source (e.g., FR I/FR II, see Czerny et al. 2009 for more 
details). 
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Table 1. The data of young radio sources. 
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G09 


44.20 


L(0 III) 


L96 


7.9 
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-1.82 


0710+439 
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kin 


P03 


45.75 


L(0 III) 


L96 


8.4 


20.0 1 


-0.75 
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1.84 
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TOO 


45.02 


L(0 III) 


L96 
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20.2 1 


-1.23 
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L(0 III) 


L96 


8.2 


19.8 1 


-1.27 
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0.007 


0.22 


kin 


L07 


45.19 


L(H/3) 


M03 


8.7 


14.6 1 


-1.69 
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1.60 
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"References for the radio sizes, the ages of these young radio sources, and the methods for deriving the ages. 
''References for the data used to derive the bolometric luminosities. 

c References for the BH masses and the apparent magnitude of the host galaxies in R— band used to estimate the BH masses in this work, of which (1) and (2) are 
selected from Snellen et al. (1996) and NED (NASA/IPAC Extragalactic Database; http://nedwww.ipac.caltech.edu) respectively. 

References. — G09: Giroletti & Polatidis 2009; V06: Vink et al. 2006; G05: Gugliucci et al. 2005; N09: Nagai et al. 2009; M99: Murgia et al. 1999; F02: Fanti 
et al. 2002; Al: this work assuming the upper limit of light speed; A93: Akujor et al. 1993; K97: Katz-Stone & Rudnick 1997; A2: this work using the jet speed 
of 0.1c in An et al. 2004; T01: Tyulbashev & Chernikov 2001; T92: Taylor et al. 1992; K06: Kunert-Bajraszewska et al. 2006; L96: Lawrence et al. 1996; M03: 
Marziani et al. 2003; X99: Xu et al. 1999; POO: Peck et al. 2000; H01: Ho & Peng 2001; X06: Xie et al. 2006; G04: Goncalves & Serote 2004; H03: Hirst et al. 
2003; W09: Wu 2009 and references therein; A00: Axon et al. 2000; W02: Woo & Urry 2002; Z03: Zakamska et al. 2003. 



